Easy- Axis Kagome Antiferromagnet: Local-Probe Study of NdsGasSiO^ 
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We report a local-probe investigation of the magnetically anisotropic kagome compound 
Nd3GasSiOi4. Our zero-field /iSR results provide a direct evidence of a fluctuating collective para- 
magnetic state down to 60 mK, supported by a wipe-out of the Ga nuclear magnetic resonance 
(NMR) signal below 25 K. At 60 mK a dynamics crossover to a much more static state is observed 
by /iSR in magnetic fields above 0.5 T. Accordingly, the NMR signal is recovered at low T above a 
threshold field, revealing a rapid temperature and field variation of the magnetic fluctuations. 

PACS numbers: 75.10.Hk, 75.30. Gw, 71.27.+a 



Kagome antiferromagnet (KAF), a 2D net of corner 
sharing triangles, is one of the most intensively stud- 
ied examples of geometrically frustrated lattices. In the 
classical case, extensive degeneracy - a hallmark of frus- 
tration - is present in the ground state (GS) manifold, 
both in the Heisenberg and the Ising limit For the 

Heisenberg case, "flexible" coplanar configurations which 
allow for zero-energy weathervane modes, are likely se- 
lected by thermal fluctuations []]]. Frustration is even 
more severe in the Ising limit, where there are minimal 
spin degrees of freedom to accommodate the compet- 
ing individual interactions. The system remains disor- 
dered at all temperatures, with exponentially decaying 
pair spin correlations [2]. In the intermediate case, var- 
ious anisotropy terms adding to the Heisenberg 
Hamiltonian can become important in selecting a partic- 
ular GS subspace. Anisotropy and additional magnetic 
field H, favoring collinear spin configurations, may lead 
to particularly rich phase diagrams, as, e.g., predicted 
for the spin-1 KAF [6j]. In this context the exciting and 
mostly unexplored field of KAF with a sizable easy-axis 
anisotropy, could well yield novel original states. 

In the recently discovered Nd 3 Ga 5 SiOi 4 (NGS) Q, 
Nd 3+ (S = 9/2) magnetic moments occupy a network of 
corner-sharing equilateral triangles, topologically equiva- 
lent to the perfect kagome lattice if the nearest-neighbor 
exchange is considered [inset (b) to Fig.[l]. Susceptibility 
curves [7] display an interesting crossover from an easy- 
plane to an easy-axis (c axis, perpendicular to kagome 
planes) behavior below 33 K, as shown in the inset to 
Fig. [2j Because of this anisotropy and large Nd 3+ mo- 
ments, NGS is unique among kagome-like compounds in- 
tensively studied in recent past [8|]. The easy-axis low-T 
behavior, originating from single-ion anisotropy and/or 
exchange anisotropy due to crystal-field (CF) effects, 
opens a route to experimental investigation of large-spin 
Ising-like KAF. The analysis of the high-T susceptibility 
yielded a Curie-Weiss temperature 9cw = —52 K, al- 
though no irregularities were observed down to 1.6 K 
Interestingly, at this low T, a magnetization plateau was 



observed at around 1/2 of the saturated (free-ion) value 
in external fields H||c of only few Tesla [7||. 

In this Letter we report a series of complementary 
local-probe measurements, including nuclear magnetic 
resonance (NMR) and muon spin relaxation (/iSR). Our 
results give a direct evidence of the fluctuating nature of 
the ground state and the absence of any static order in 
zero magnetic field. In addition, we investigate a dynam- 
ics crossover of the magnetic fluctuations when H||c is 
applied, toward a more static polarized state. Our find- 
ings prove persistent spin fluctuations in the zero-field 
ground state down to 60 mK, characteristic of a spin- 
liquid state, which has been recently suggested from the 
absence of neutron magnetic Bragg peaks and the pres- 



2 







jL"c(3) W 


3? 






Ga(3)/Sifv 




boo) 


Int. (a 










fw v> 


§ 
Z 






S "C(3) 






3.6 3.7 3.8 
H(T) 




J\ 








j 






1.2 K 








'. t \— 

If W_ 

? S< P 




— «- — f f t — . i . — 

, i l S(3)™S(3) ,1 S(l) "S(1)"'S(1)"S(3) 


"sp> i 


! 3 






4 


5 


6 



H(T) 



FIG. 1: (color online). 69,71 Ga NMR spectra of NGS at 
46.01 MHz and H||c. The 1.2 K spectrum is translated verti- 
cally. Arrows indicate the satellite (S) lines at 80 K. The inset 
(a) highlights central (C) peaks of the three crystallograph- 
ically non-equivalent Ga sites at 80 K; the vertical dashed 
line corresponds to a reference field. The inset (b) shows 
the network of Nd 3+ magnetic moments and the location of 
the in-plane Ga(l) sites and Ga(2, 3) sites separating kagome 
planes. 
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FIG. 2: (color online). Scaling of the 71 Ga NMR shift with 
bulk susceptibility. Inset: The susceptibility crossover at 1 T. 



ence of a weak diffuse scattering down to 46 mK jg] . Par- 
tial ordering has been also reported under applied field 
at low T. 

Polycrystalline NdsGasSiO^ samples were synthesized 
by a solid state reaction and a single crystal of cylin- 
drical shape (2 mm in diameter and 8 mm in length) 
was grown by a floating zone technique 0. NMR mea- 
surements were performed at various fixed frequencies by 
sweeping the magnetic field. Nuclear quadrupolar reso- 
nance (NQR) was also measured at zero field. Powder 
/iSR measurements were performed on the LTF and GPS 
spectrometers at PSI, Switzerland. 

We first report our Ga NMR investigation. The lo- 
cal magnetic coupling with the Nd 3+ moments is dif- 
ferent for the three crystallographically non-equivalent 
gallium sites Ga(l, 2, 3) [inset (b) to Fig. Q]. In addi- 
tion, both 69 ' 71 Ga isotopes contribute to the NMR spec- 
trum. Through their sizable quadrupolar moments they 
are sensitive to local charge distributions, creating crystal 
electrostatic field gradients (EFG). We assigned all the 
observed peaks shown in Fig.Q]- three sets of central (C) 
and two satellite (S) lines for each isotope - by perform- 
ing frequency and angular dependent measurements. We 
took advantage of the EFG tensors symmetry [13], known 
gyromagnetic ratios ( 69 j/ 71 j — 0.7871) and the ratio of 
the quadrupolar moments ( 69 Q/ 71 Q = 1.583). Since the 
Ga(3) site is randomly occupied by Ga 3+ and Si 4+ 0], a 
large distribution of quadrupolar frequencies is expected. 
This results in broad S lines and a broad NQR spectrum 
ficj |- Surprisingly, the distribution of local environments 
does not lead to glassy-like behavior at low T . 

In order to eliminate the quadrupolar broadening and 
get the best accuracy on the local magnetic proper- 
ties, the magnetic field was applied perpendicular to 
the kagome planes (H||c), along a local three-fold ro- 
tational axis for the Ga(l, 3) sites. This results in 
narrow C(l, 3) lines, on which we now focus. Their 



shift K = (Hq — H c )/H c = A\i + Kq of a resonance 
field H c from a reference field Hq is proportional to 
the local susceptibility xi- It has three contributions, 
which need to be considered in the coupling constant 
A = Ahf + A dd + A dm ; a hyperfine (hf) and a dipo- 
lar (dd) coupling between Nd 3+ moments and Ga nuclei, 
and a macroscopic demagnetization (dm) field. Kq is 
a T-independent chemical shift. Our calculation of the 
dipolar field [H| yields that all these contributions are of 
a similar magnitude and that the hyperfine coupling is 
almost isotropic [l(J. In Fig. [2] we show the linear scaling 
of the shift with bulk susceptibility between 80 K and 
25 K. This proves that the local and the bulk magnetic 
response are identical and, in particular, the 33 K cross- 
ing is also detected by our local-probe measurement. On 
the contrary, in formerly studied 3d metal kagome-like 
compounds the intrinsic susceptibility is regularly over- 
shadowed by magnetic defects contributions at low T [§]. 

Next, we investigate the Nd 3+ spin dynamics. In 
Fig. [3ta) we show the T-dependence of the 71 C (3) 
spin-lattice relaxation rate 1/Tj for H||c. Its val- 
ues were obtained by fitting a magnetization recov- 
ery curve after saturation to the relaxation func- 
tion M z (t) = M z [1 - 0.4exp(-t/Ti) - 0.6exp(-6t/7i)], 
suited for magnetic relaxation of the central line in spin- 
3/2 nuclear systems [l3|. The relaxation is rather fast, 
proving that it is due to rapidly fluctuating dipolar and 
hyperfine fields, originating from Nd 3+ moments. With 
lowering T the magnetic fluctuations are found to slow 
down considerably, which significantly enhances 1 jT\ be- 
low 100 K and wipes out completely the signal below 
25 K. The wipe-out is depicted in Fig [31(b), displaying 
the T-dependence of the NMR intensity. It is caused by 
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FIG. 3: (color online), (a) Zero-field (ZF) muon spin re- 
laxation rate A and 71 C (3) NMR spin-lattice relaxation rate 
1/Ti measured at 46.01 MHz (3.55 T) and H||c. Inset: ZF 
relaxation of muon polarization, (b) Wipe-out of the NMR 
intensity. The high-T (low-T) data were taken at 46.01 MHz 
(29.61 MHz) at different fields. Dashed lines are a guide to the 
eyes. Inset: Absence of the NQR echo at 1.2 K and 15 MHz. 
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the increase of both the spin-lattice and the spin-spin 
relaxation rate. In a conventional transition scenario to 
a frozen state one would expect the NMR/NQR signal 
to emerge again below the transition temperature. In 
contrast, the inset to Fig.[3](b) demonstrates the absence 
of the NQR signal at 1.2 K, which points to persistent 
magnetic fluctuations. Similarly, the NMR signal is not 
detected at 1.2 K in small applied fields (H < 3 T). 
Only in higher fields the NMR signal is recovered [Fig. [IJ 
Fig-fH^b)]. We will address this field-induced effect later. 

Further insight into the Nd 3+ dynamics is provided by 
zero-field (ZF) fiSR measurements. The muon relaxation 
rate A could be followed in the full T-range from 100 K 
down to 60 mK and was extracted from a stretched- 
exponential fit P(t) = Poexp [— (Xt) a ] of the muon po- 
larization. The T-independent a = 0.6(1) likely reflects 
a variety of oxygen, hence fi + sites in NGS. In the sim- 
plified case of an exponential relaxation, A is related to 
the Nd 3+ magnetic fluctuation rate v e by [13] 
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FIG. 4: (color online). Field dependence of the muon relax- 
ation rate A at 60 mK and the 69 Ga NMR 1/Ti at 1.2 K 
on several C and S lines. The solid line gives the predicted 
-H-dependence of A in the dynamical case [Eq. fT}]. Inset: Re- 
laxation of muon polarization in various LF fields at 60 mK. 
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where H^ = 0.2 T is the average dipolar field created by 
the Nd 3+ spins at the muon sites, 7 M /2?r = 135.5 MHz/T 
is the muon gyromagnetic ratio and Hlf stands for the 
applied longitudinal field if any. The thermal evolution 
of A oc l/v e in zero field, plotted in Fig. [3](a) reflects 
a strong slowing down of the Nd 3+ spin fluctuations - 
from ~ 200 GHz at 100 K, in agreement with paramag- 
netic fluctuations governed by the magnitude of the ex- 
change, to ~ 4 GHz at 10 K. The saturation of A below 
10 K indicates that the slow magnetic fluctuations, likely 
associated to strongly short-range coupled Nd 3+ spins, 
persist in NGS down to 60 mK, which is a fingerprint of 
a d yna mical ground state, observed in many kagome- like 
0,13,0, EE El and pyrochlore (H ED, El lattices. In 
most of these materials the dynamical plateau though co- 
exists with a spin-glass- like transition or even LRO. Such 
transitions are regularl y ob served at the temperature of 
the plateau onset [H, E3, E3, El, ED, HI- Supported 
by a lack of anomalies in the magnetic susceptibility of 
NGS 0, the muon relaxation curves, which relax to zero 
down to 60 mK [inset to Fig. [3|a)], provide a direct evi- 
dence of the absence of any static magnetic fields emerg- 
ing from the Nd 3+ moments. Our /iSR data thus sup- 
port a collective paramagnetic state in NGS, stabilized 
below 10 K, which makes it together with the recently 
discovered spin- 1/2 Heisenberg system ZnCu3(OH) 6 Cl2 
[23l | the only KAF potentially displaying a pure spin- 
liquid ground state [24]]. 

When applying the longitudinal field (LF) at 60 mK, 
A varies moderately up to 0.5 T and then suddenly 
decreases by more than two orders of magnitude for 
Hlf = 2.5 T, as visible already in raw relaxation curves 
(Fig. [4]). Such a strong response to Hlf, not observed 



in any other kagome- like compound 16|, 17, 11, cannot 
be explained within a fast fluctuation regime. Eq. (Q]) 
predicts a smooth and much slower decrease of A with 
the upward shift of the Larmor frequency u>l = J^Hlf 
with respect to the fixed fluctuation rate u e [solid line in 
Fig. [4]. The observed decrease of A above Hlf > 0.5 T 
demands that the spin dynamics itself, v e , is strongly 
reduced by Hlf, which suppresses the spectral density 
of fluctuations at u>l- Assuming if M = 0.2 T, Eq. (U) 
gives v e ~ 10 MHz at 2.5 T. Therefore, our data point 
to a field-induced crossover from a dynamical to a more 
static state for Hlf > 0.5 T at 60 mK, in agreement 
with the recent detection of a static, net-ferromagnetic, 
state in the same T and H range d- 

The field-suppressed magnetic fluctuations are further 
reflected in a pronounced JT-dependence of the NMR 
spin-lattice relaxation rate 1/Ti at 1.2 K, as measured 
at various fields in the spectra and at various frequen- 
cies (Fig. |4]). 1/Ti is lowered by almost three orders 
of magnitude in a narrow JT-window between 3 T and 
6.5 T. This enhanced effect of H on the magnetic fluctu- 
ations is responsible for only partial observability of the 
NMR spectrum at low T. The 3.5 T threshold field at 
1.2 K and 46.01 MHz is easily extracted from the miss- 
ing S(l, 3) lines below this field, as evidenced in Fig. [TJ 
The observed peaks are significantly shifted and broad- 
ened with respect to the high-T spectrum. Using the 
shift of the narrowest 71 C(1) line and its high-T coupling 
constant A, the Nd 3+ magnetic moment of 1.6 [is, ca. 
half a value of the full moment, is extracted. Note that 
this is the mean value of the six moments on the kagome 
star, to which the Ga(l) is coupled. Individual Nd 3+ mo- 
ment could be larger than 1.6 /is in a non-ferromagnetic 
arrangement on the star. This value is constant above 
the threshold field, in accordance with previous magne- 
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tization measurements Q. The static line-shift and the 
dynamical relaxation thus exhibit strikingly different in- 
dependence at 1.2 K. When lowering the temperature be- 
low the wipe-out the NMR signal reappears in a rather 
narrow low-T range. In Fig[3](b) we show that the inten- 
sity is completely recovered between 1.7 K and 1.2 K at 

4.5 T, while at lower fields it is not fully recovered yet 
at 1.2 K. This points to a rapid temperature variation 
of the threshold field, which is a strong indication of a 
field-induced transition in NGS. 

On the theoretical side, considerable attention has 
been recently devoted to the magnetic-field effects on the 
ground state of the Ising-like KAF, either due to easy- 
axis single-ion anisotropy (SIKAF) 0, @| or anisotropic 
exchange (AEKAF) [2 HE 0|- Since the CF sym- 
metry on the Nd 3+ site is low in NGS Q, one ex- 
pects that the lowest Kramer's doublet essentially de- 
scribes the low-T physics. A priori, one can not dis- 
card single-ion nor exchange anisotropy. If the exchange 
J = 3\0cw\/4:S(S + 1) w 1.5 K is much lower than 
the CF splitting between the ground and excited states 
the SIKAF picture will apply, otherwise the exchange 
anisotropy may be important. For this latter AEKAF 
case, an "exotic ferromagnetic" phase, showing finite 
magnetization but no conventional order, was predicted 
at H = with a transition temperature T c < 0.08J Z S/A, 
which depends on the asymmetry level A = J z /J X y 0|- 
To comply with the absence of transition above 60 mK in 
NGS, this model would yield a strong easy axis scenario 
with A > 10 which does not fit well with the moder- 
ate magnetocrystalline anisotropy observed at low T 0. 
NGS is thus likely best described in the SIKAF model 
with finite single ion anisotropy D. The susceptibility 
crossing at 33 K suggests that D is of the order of 10 K. 
Our NMR results at 1.2 K together with magnetization 
data at 1.6 K 0] point to a magnetization plateau with 

1.6 [ib per Nd 3+ . It is worth noting that the expectation 
value of the total angular momentum operator within 
the ground Kramer's doublet is likely reduced with re- 
spect to the full free-ion value (3.27 fin)- Therefore, the 
magnetization fraction M/Mf on the observed plateau 
is not accurately known, but M/Mf > 1/2. In striking 
contrast, all anisotropy models on the kagome lattice pre- 
dict a 1/3-magnetization plateau [^HEHol- This plateau 
should be stable up to a field of ca. 10 T (H ~ JS) 0, 
in contrast with the experiment. 

In conclusion, our results point to the persistence of 
the fluctuating disordered state in NGS down to 60 mK 
in zero field. The rather slow fluctuations, v e ss 4 GHz, 
witness a collective paramagnetic state below 10 K, which 
is a new relevant energy scale in the system. NGS and 
ZnCu3(OH) 6 Cl2 thus enable an experimental verification 
of theoretical predictions of a spin-liquid ground state 
for quite different Hamiltonians. The fluctuations in 



NGS are considerably suppressed by the applied field 
and signify an unpredicted field-induced phase transition 
above 0.5 T at 60 mK. Much like in rare-earth based py- 
rochlores, our results speak in favor of competing single 
ion anisotropy and exchange interaction. Such a compe- 
tition may lead to the observed collective paramagnetism 
and to its instability against the applied magnetic field. 
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